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We present a new dispersive analysis of the isospin breaking decay r\ 37r. The resuhing rep- 
resentation of the decay amplitude allows us to determine the quark mass double ratio Q and we 
find as a preliminary result Q = 22.3 ± 0.4. Finally, we discuss a number of improvements that 
we intend to implement in the future. 
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1. Introduction 

The decay rj — > 371 is forbidden by isospin symmetry, as Bose statistics does not allow three 
pions to form a state with vanishing total isospin and total angular momentum. Neglecting elec- 
tromagnetic contributions that are strongly suppressed [1], the decay amplitude A(5,f,M) is propor- 
tional to {nid — niii) or, alternatively, to the quark mass double ratio 

1 m^ — m^, 

— = — (\ \) 

^ — 

where m = {rriu + md) /I. As F oc |Ap oc Q we can get Q by comparing a theoretical result for 
the amplitude with a measurement of the decay width T. Alternatively, Q can also be calculated 
from a ratio of meson masses [2, 3]. 

As is well known, the chiral perturbation theory series converges rather slowly: At tree-level, 
the decay width is F = 66 eV [4, 5], while at one-loop it is already F = 160 eV [6], which is still 
quite far away from the experimental value F = 295 eV [7]. Of course, the theoretical values need 
an input for {m^ — niu), which is calculated from meson masses: 

{nid - mu)Bo = (m^o - )ecD = 'w^o - m\+ - m^o + m]^+ , (1.2) 

where the second equality relies on Dashen's theorem, stating that at leading order in the low energy 
expansion the electromagnetic contribution to the mass differences m^o — mj^+ and — m^o is 
the same [8]. The slow convergence of the chiral series is mainly due to final state rescattering of 
the pions, which can be very well treated by means of dispersion relations. 

The idea to use dispersion techniques to calculate the rj ^ 37i amplitude is, in fact, not new: 
it has been done akeady more than 10 years ago by Kambor, Wiesendanger and Wyler [9] and 
by Anisovich and Leutwyler [10]. The methods used in these works differ in technical aspects, 
but lead to compatible results. We follow the approach as presented in ref. [10]. There has been 
a lot of activity in this area since then and therefore, we think it is worth to take a fresh look at 
this problem. There has been considerable improvement concerning the nn phase shifts [1 1 - 14], 
which are the most important physical input in the dispersion relations. Furthermore, there are a 
number of new measurements of this decay by the KLOE [15, 16], MAMI [17, 18] and WASA 
[19, 20] collaborations, which will be useful for the determination of the subtraction polynomials, 
and also a full two-loop calculation in chiral perturbation theory [21]. 



2. The Method 

We only give a very brief account of the method that we are using and refer to ref. [10] for 
a more detailed description. For the charged decay rj — > tt^tt+tt^, we define the Mandelstam 
variables as 5 = {pj,:+ + p^-Y, t = {pn- + P^ro)^ and u = {pj^ + Pn+Y- They are related hy s + t + 
u = m^ + m^o + 2m^+ = 3^0. At the order we are working the amplitudes for the charged and the 
neutral channel aie related by 

^neutral('^'^'") ~'^charged('^'^'") ~'"'^charged(^'"''^) ~'"'^charged("''^'0 • (2-1) 
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In the following we will restrict ourselves to the charged channel, as we could get the neutral 
channel easily by the above relation. We are working in the isospin limit where m^o = m;j+ . It is 
not a priori clear, what value for the pion mass in the isospin limit should be used and we choose 
= nij^Q. 

Unitarity allows us to relate the imaginary part of the decay amplitude An=A^^n, where the 
final state n is some three pion state, to the nn scattering amphtude and A„ itself as 

ImA„ = {2nf5\pn-Pn')T:,, }A„, . (2.2) 

n' 

This is a linear constraint for the amplitude and we can thus extract a normalisation factor 

A{s,t,u) = --^ J' M{s,t,u). (2.3) 

With this choice, the tree-level expression for M{s,t,u) is normalised to 1 at the centre of the Dalitz 
plot where s = t = u = sq, and the one-loop amplitude only depends on measurable quantities. In 
particular it does not depend on Q. 

If the discontinuities of D- and higher waves are neglected, we can decompose the amplitude 
into isospin components in the same way as was done for the Tin scattering amplitude in ref. [22]: 

M{s,t,u) =MQ{s) + {s-u)Mi{t) + {s-t)Mi{u)+M2it)+M2iu)-^M2{s). (2.4) 

As only the S- and P-waves have discontinuities up to two-loop order in chiral perturbation theory, 
the decomposition is exact up to that order. 

Inserting this decomposition in eq. (2.2), we get for the discontinuities of the isospin ampli- 
tudes 

disc Mj{s) = Mi{s + ie)-Mj{s-ie) ^ ^^^^^^ +M;(.)} .-'^'^ sin5,(.), (2.5) 

where 5j{s) are the S- and P-wave nn scattering phase shifts and / = 0, 1,2. The inhomogeneities 
M]{s) consist of angular averages over the Mj: 

2 20 2 

Mo{s) = -{Mo)+2{s-so){Mi) + —{M2) + -K{m), (2.6a) 

Mi{s) = i |3(zMo) + ^{s-so){zMi) - 5{z.M2) + ^KizM)^ (2.6b) 
M2(^) = (Mo) -^{s- so) (Ml ) + ^ (M2) - ^ K{zMi ) , (2.6c) 



where 



and 



s — 4ml. 



s 



k{s) = \ ^ \ {{mn ^mnY - s]{{mr^ - m^Y - s] . (2.8) 



The integration variable in eq. (2.7) is z = cos 0, where Q is the scattering angle. 
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These expressions for the discontinuities enable us to write down a set of dispersion integrals, 
coming from Cauchy representations of the functions Mj{s)/D.j{s) [10]: 



Mo{s) = D.o{s) < 



. f ds' sin5o(/)Mo(/) I 



4ml 



\ao{s')\{s' - s - ie) ( ' 



/ ^ ^ , . \ s f ds' sin 5i (s' )M I (s') I 

4m2 



/ ^ ^ , s t ds' sm&i{s')M2{s') 

Q.i[s) are the so called Omnes functions [23], which are the solutions of eq. (2.5) for Mi{s) = 
and are given by 



Q.i{s) = exp < 



' ' ds'y (2.10) 



K J s'{s' — s 

4ml 



Eqs. (2.9) contain four parameters Oq, jSq, Jo and jSi. These subtraction constants are not de- 
termined by the dispersion relations and have to be fixed otherwise. Actually, there are three more 
of these constants in the equations for Mi and M2, but they can be eliminated because the decom- 
position given in eq. (2.4) is not unique. The Mj can be shifted by a polynomial without changing 
the total amplitude M{s,t,u), and the polynomial can be chosen such that three subtraction con- 
stants vanish. The remaining four are determined by a matching to the one-loop result from chiral 
perturbation theory. 

We solve the integral equations numerically by an iterative procedure: We set the Mj{s) to the 
tree-level result from chiral perturbation theory and calculate the Mi{s) by eqs. (2.6) and then get a 
new result for the Mj{s) from eqs. (2.9). The subtraction constants Oq and jSq can now be calculated 
by a matching to the one-loop amplitude at the Adler zero, as the amplitude there is protected by 
the chiral SU(2) x SU(2) symmetry and does not change much after the one-loop level. While 
these two constants have to be determined in every iteration step, the other two, 70 and jSi, can be 
set to their final value right from the start. Once the four constants are fixed, we add the subtraction 
polynomial and then repeat the same procedure until the subtraction constants Oq and /3o converge. 

Special cai^e has to be taken when evaluating the angular- averages M]{s) in eqs. (2.6), as the 
starting and end points of the integration path, which are functions of s, can take complex values. In 
the first iteration step, the functions M[{s) are analytic everywhere and the integral does not depend 
on the path. However, the situation changes afterwards: the Mi{s) have a cut along the real axis 
and the integration path has to be deformed in such a way that it does not cross the cut [24] . 

3. Preliminary Results 

In the future, we will incorporate a number of improvements to the method described in the 
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Figure 1: The function M(){s) and its change during the iterations. The left panel shows the real part, the 
right one the imaginary part. The red dotted line is the tree-level ;i^PT result used as initial configuration, 
the blue dashed and the green dash-dotted lines are the results after one and two iterations, respectively. The 
solid black line is the final result after convergence of the subtraction constants is reached. In the real part, it 
is almost impossible to distinguish the result after two iterations from the final result. 
To have a smooth representation of Mo{s) near s — 4m^, we have to ensure that the threshold of the phase 
shifts agrees with the lower limit of integration. As we know the phase shifts only for m„ = we also 
used the charged pion mass to create these figures. 



last section. The results we present here are therefore only preUminary. The plaimed extensions 
are discussed in the next section. 

Fig. 1 demonstrates how the function Mq (5) develops during the iteration procedure. While for 
the real part the curve does not change considerably after the second iteration step, the imaginary 
part shows that it is worthwhile going beyond that, as there is still an obvious difference between 
the result after two iteration steps and the final result. 

From the functions M/(s) we can then calculate the amplitude and the Dalitz plot. We use the 
dimensionless standard Dalitz plot variables defined as 

where 

Qri = '"77 - 2m^+ - m^a . (3.2) 

Fig. 2 shows our result for the Dalitz plot, normalised to 1 for X = 7 = 0. It agrees with the 
polynomial representation provided by KLOE [15] at the 10 % level. 

With the experimental decay width F = 295 ±20 eV [7], we get for the quark mass double 
ratio Q = 22.3 ± 0.4. This is in good agreement with other results presented in fig. 3, which range 
from 20.7 to 24.3. The error is only due to the experimental uncertainty in the decay width, as we 
do not have an estimate for the theoretical error yet. 

As a last illustration, fig. 4 shows the development of Q over several iteration steps. Already 
after two iteration steps, the result for Q is well within the experimental error band of the final 
result. However, going beyond that will not only shift the central value within the error bounds, but 
also reduce the theoretical error, which is not included in this figure. 
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Figure 2: The Dalitz plot for tj n^n^n normalised to 1 for X = 7 = 0. The dimensionless DaUtz plot 
variables X and Y are defined in the text. 
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Figure 3: A selection of results for Q. Our result is 22.3 ± 0.4 and is indicated by the grey band. The error 
is only due to the experimental uncertainty on the decay width. The other results are taken from Leutwyler's 
talk at this conference [25], from a dispersive analysis in ref. [9], from a two-loop calculation in xPT^ [21], 
from Weinberg's quark mass ratios [2] and from an analysis including Dashen violation [3]. The last value 
we calculated from the MILC quark mass ratios presented by Heller at this conference [26]. 
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Figure 4: The development of Q over several iteration steps. The blue line is the final result after many 
iterations, the blue band is the corresponding error coming only from the uncertainty of the decay width. 



4. Conclusion and Outlook 

We have presented a new dispersive analysis of the decay rj 3n and argued that there have 
been several important developments in this and related fields that make a new analysis worthwhile. 
A number of preliminary results have been given, in particular we get for the quark mass double 
ratio Q = 22.3 ± 0.4, which agrees well with a number of results from other works. We stress again 
that this number is preliminary and that the error given is only due to the uncertainty in the decay 
width and does not contain any estimate for the theoretical error. 

The method as presented here has a number of shortcomings that we will have to address in 
the future. Besides solving these problems, we will also include additional features that have not 
been part of earlier dispersive treatments of this decay and that the precision reached in experiments 
nowadays has rendered mandatory. 

As mentioned above, we calculate the decay amplitude in the isospin limit and use the mass 
of the neutral pion. However, it is not a priori clear, what value for the pion mass should be used 
and indeed the final result depends on this choice: setting the pion mass to mj^+ shifts Q from 22.3 
down to 21.0. The pion mass enters the dispersion integrals in eq. (2.9) directly, as it is contained in 
the integration limits and in k{s) , and also via the phase shifts. While it is easy to use the physical 
masses for the former, it is not for the latter. We would need different phase shifts for scatterings of 
pions with different flavours, which are not available at the moment. Indeed the above-mentioned 
shift in Q is obtained without altering the pion mass in the phase shifts. 

A number of other contributions have not been included yet: electromagnetic interactions 
[27, 28], inelasticity and the imaginary parts of D- and higher partial waves. Future extensions of 
the algorithm will include estimates of these effects. 

The subtraction constants have been determined so far by a matching to the one-loop result 
from ;f PT. We plan to make use of experimental data instead and calculate the subtraction constants 
by fitting our Dalitz plot to a measurement. 
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Of course we can and will also extract a number of other interesting parameters than just Q, e.g. 
the quark mass ratio R, the branching ratio r = r^j,.o jT^fi^+j^- or the quadratic slope parameter a. 
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